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ABSTRACT 

Abundance of rich clusters in local universe is currently believed to provide the most ro- 
bust normalization of power spectrum at a scale of 10 Mpc. This normalization depends 
very sensitively on the calibration between virial mass M and temperature T, which is 
usually taken from simulations. Uncertainties in the modelling, such as gas cooling and 
heating, can lead to a factor of two variations in the normalization and are thus not 
very reliable. Here we use instead an empirical M500 — T relation derived from X-ray 
mass determinations to calibrate the method. We use results from dark matter simula- 
tions to relate the virial mass function to the mass function at observed M500 • We find 
that the best fitted value in flat models is erg = (0.7±0.06)(O^/0.35)-°-44(r/0.2)0 08, 
where only the statistical error is quoted. This is significantly lower than previously 
obtained values from the local cluster abundance. This lower value for erg is in a bet- 
ter agreement with cosmic microwave background and large scale structure constraints 
and helps alleviate small scale problems of CDM models. Presently the systematic un- 
certainties in the mass determination are still large, but ultimately this method should 
provide a more reliable way to normalize the M — T relation. This can be achieved by 
obtaining a larger sample of well measured cluster masses out to a significant fraction 
of virial radius with BeppoSAX, Chandra and XMM-Newton. 

Key words: cosmology: theory - dark matter - galaxies: haloes - galaxies: clusters: 
general - X-rays: galaxies. 



1 INTRODUCTION 

Abundance of clusters in the local universe provides one 
of the strongest constraints on the matter power spectrum 
amplitude on scales of order lOMpc. The main power of this 



data sets also provide some more model dependent con- 
straints on the same scale. Recent determinations of the best 
f itted model from CMB combined wi t h large scale structur e 
( |Wang, Tegmark, fc Zaldarriaga 200l|, [Efstathiou et al. 200ll) 



methoci is LliuL Lliy abuudaucc depyuds vyiy syuslLlvyly un 
the am plitude of linear power apcetrum on thia acalc. I'Ju - 
merous studies over the past decade have investigat ed this 



and from Ly-a forest (Croft et al. 199!:, [McDonald et al 



Liddle 1996], [Kitayama fc Suto 19971 |Eke et al. 199^, [P~ 



constraint, obtaining more or less consistent resul ts ( White 
Efstathiou, fc Frenk 1993t [Bond fc Myers 1996|, |Viana & 



200c ) . Two mos t recent studies are those by Pierpaoli, Scott 
fc White (2001 ) using a local sample of clusters and [Borganl 



199a , [Wang fc Steinhardt 1998| , [Blanchard et al. 2000|, |Hcn~ 



200c) both give ag, values more than 30% lower than those 
from the latest local cluster abundance determinations for 
the favored fim = 0.3 — 0.4 case. A lower normalization also 
helps si gnificantly with the problems that CDM has on small 
scales (Alam, Bullock, fc Weinberg 2001). This tension be- 



et al. (2001) using a higher redshift sample, both of which 
find a relative error of 5 — 10% on the density dependent 
normalization (but which are inconsistent with each other 
at a more than 3-(t level, see discussion below). 

While cluster abundance normalization is usually as- 
sumed to be the most reliable on a scale of lOMpc, other 
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tween the cluster abundance normalization and other data 
sets, which is formally at more than 3-(t level, could indicate 
there are still some unrecognized systematic uncertainties 
present in the cluster abundance normalization. 

An implicit assumption in all of the studies above is to 
assume a relation between halo mass M and observed gas 
temperature T, M / {lO^^h-^ Mq) = (T//?)^/^ A^^''^S-\ 
where Ac is the mean overdensity inside the virial radius 
defined using spherical collapse model, T is gas temperature 
in keV, — Q,m{)-+ z)'^ +0.^ for the zero-curvature case and 
the constant /3 is determined from the hydrodynamical sim- 
ulations. The problem with this approach is that the simula- 
tions may not capture all of the physics occuring in the real 
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world. For example, most simulations only include adiabatic 
hydrodynamics, while it was recently shown that including 
gas cooling (which must be present at some level) can lead to 
a factor of two change in th e normalization constan t of the 
mass temperature relation ( Muanwong et al. 2001 ). Other 
potentially important effects, such as heating of the gas with 
stars, supernovae or AGN, could further change the relation. 
Resolution issues could also be affecting the results: there is 
quite a lot of scatter in the relation between different gr oups 
even w hen all the phy sical ingredients are ident ical (see 



Ko- 



matsu .'z Seljak 2001 and Pierpaoli et al. 2001 for a recent 
comparison between different groups). Yet another possi- 
ble complication is the observer's definition of temperature, 
which depends on the spectral bandpass a nd other instru- 
mental details. Mathiesen & Evrard (2001) explore several 



different possibilities and find up to a 20% variation between 
them. Since all these effects are systematic it is difficult to 
place a reliable error estimate on them, except to argue that 
a factor of two changes are not excluded at this point. To 
put this into a context, such a change leads to a 30-40% 
change in the overall power spectrum normalization, much 
larger than the quoted error-bars. It seems therefore worth- 
while to explore alternative methods to normalize the M — T 
relation which do not rely exclusively on the simulations. 

Over the past several years there has been a lot of 
progress in the direct determination of cluster masses using 
spatially resolved temperature and X-ray intensity informa- 
tion. Under the assumption of spherical symmetry and hy- 
drostatic equilibrium with no non-thermal pressure support 
one can solve for the radial mass profile of the matter in the 
cluster. This approach consistently yields M — T relation up 
t o a factor of two lower than the ave r age over simulations 
( Horner, Mushotzky, fc Scharf 199i^ , iNcvalaincn, Marke- ' 



vitch, Forman 2000, Finogucnov, Reiprich, fc Bohringci 



2001, Mien, Schmidt, fc Fabian 2001). Such a discrepancy 



could be explained if there was an additional non-thermal 
pressure support, which should be added to the equation of 
hydrostatic equlibrium (such as from turbulent motions or 
magnetic fields). However, recent comparison between X-ray 
mass measurements and weak lensing mass measurements 
for several clusters shows very good agreement between the 
two, excluding the possibility of a large non-thermal sup- 
port, at least for relaxed clusters (?Vllen et al. 2001). An 



alternative possibility is to assume that X-ray mass deter- 
minations are reliable and there is some systematic problem 
with the simulations. 

The purpose of this paper is to use the empirical M — T 
relation to derive the power spectrum normalization. An ear- 
lier attempt to do this based on one m easured cluster mass 
has been made by Markevitch (1998) and a lower value of 



erg has been found than calibrating from simulations. Since 
then the observed M — T relation and scatter around it has 
been established much more accurately. Since the measured 
masses are limited to the inner parts of the cluster they can- 
not be extended to the virial mass, defined here as the mass 
that gives rise to the universal mass function, without ad- 
ditional modelling. Here we use CDM type mass profiles to 
make this connection. We find that the offset relative to the 
simulations persists and leads to a significantly lower power 
spectrum normalization than was found before. The impli- 
cations of this result and future prospects for this approach 
are discussed in the conclusions. 



2 THEORY 

The halo mass function describes the number density of ha- 
los as a function of mass. It can be written as 

dn _ P r( . dlna^^ 
dlnM ^ hr^'"' d\nM ' ^ ' 

where p is the mean matter density of the universe, M is 
the virial mass of the halo and n{M) is the spatial number 
density of halos of a given mass M. We introduced func- 
tion /((t), which has a universal form independent of the 
power spectrum, matter density, normalization or redshift if 
written as a function of rms variance of linear density field 



a{M)=^-K / P{k)WR{k)k^dk. 



(2) 



Here WR{k) is the Fourier transform of the spherical top hat 
window with radius R, chosen such that it encloses the mass 
M = A-kR^p/Z and P{k) is the linear power spectrum. 

The universality of the mass fun ction has been recently 



investigated by a number of auth ors (Sheth & Tormen 199!: , 
Jenkins et al. 2001, White 2001), where it has been shown 



that the mass function is indeed universal for a broad range 
of cosmological models. Jenkins et al. (2001) propose the 
form 



/(cr) = 0.315 exp[- 



Ina"' +0.61|^-*1 



(3) 



They find this equation works best if the cluster mass M 
is defined as the mass within the radius where m ean den- 
sity in units of critical is 200S7m. Other authors (Sheth & 



Tormen 1999| , [White 200 1| ) find that a similar mass function 
works if the mass is defined within the spherical overden- 
sity defined by the spherical collapse model, which is 108 
for ACDM model with fim = 0.35 th at we adopt here a s 
the reference model. We will adopt the Jenkins et al. (2001) 
value here and comment on the implications of other choice 
in the discussion section. 

Given that the mass profiles at large radii are difficult to 
measure observers prefer to quote the masses Ma at smaller 
radii or higher A than viri al values. Two of ten used values 
in the literature are Af25oo (Allen et al. 2001) and M500 (e.g. 



Finoguenov et al. 2001). Since at the moment we do not have 
mass functions at M500 or M2500 directly from simulations 
we will construct them by relating Ma to Afvir adopting 
an average mass profile for the clusters. This is typically 
parametrized as (Navarro, Frenk, & White 1997) 



p{r) 



(r/rs)-«(l + r-/rs) 



3-l-a ■ 



(4) 



This model assumes that the profile shape is universal in 
units of scale radius Ts, while its characteristic density ps 
at Ts or concentration ca^ = r Mali's rnay depend on the 
halo mass. The halo profile is assumed to scale as in the 
outer parts and as r" in the inner parts, with the transition 
between the two at r s. We fix the inner slope to a = — 1 
(Navarro et al. 1997), since we are not concerned about 



the shape of the halo profile in the center. The outer slope 
is the most common value found in the A'^-body sim- 
ulations, al though scatter arou nd this value can be quite 
significant (Thomas et al. 2001). Even if we adopt this value 
the profile is still not fixed and the remaining freedom can 
be parametrized with the concentration parameter ca^. In 
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would find (3 = 1.95 ± 0.25. This value should be com pared 
to /3 = 1.3 =b 0.2 as the average over the simulations (Pier- 



paoli et al. 2001). The two values are inconsistent with each 



other at a 3-cr level, suggesting that there is a significant 
systematic error either in the observed mass estimates or in 
the simulations. Other data samples find similar discrepan- 
cies between simul ated and observed relation. For example, 
Allen et al. (200l| ) find a 40% offset in M2500 - Tbsoo re- 
lation relative to the simulations. It is unclear how r2500 
obtained from Chandra relates to temperatures in cluster 
samples used to derive the power spectrum normalization, 
which are mostly obtained from ASCA. Moreover, M2500 to 

conversion carries considerable uncertainty (figure 
For these reasons we will not use this data set here. 

For the cluster s ample we use a recent compilation by 
Picrpaoli et al. (2001) based on the data from nearby clus- 
ters. Their adopted tem perature measurements are taken 
mostly from White (2OOOI), which are on average 15% higher 
than those derived by Markevitch (199q ). This by itself can 
be a significant systematic effect. Since w c wish to compare 



this d ata sample to the M — T relation by Finoguenov et al 



(2001) we must compare the temperatures of the two sam- 



Figure 1. M/^/M^^ versus A for = 3.7 (short dashed), 

CAa = 5 (solid) and ca^. = 7.5 (dotted) for Ac = 108. Also shown 
is isothermal profile Ma/Ma^ = (A/Ac)^^/^ as long dashed 
line. 



pies. We do this in two different ways. First, we compare the 
derived luminosity-temperature (L — T) relation for the two 
samples, which gives us the relative normalization between 
the adopted temperature values, si nce in both samples lumi - 
nosities are obtained from ROSAT. Finoguenov et al. (2001 ) 
data combined with L ~ M relation (Roiprich & Bohringer 



sity in units ot AiAc tor ca^ = d.Y,t), Y.O, which spans the 
range of concentration parameters appropriate for massive 



figure g we show mass as a function ot spherical overden- gi^es L = 1.15 x 10'*'^(r/6keV) ^fe~^ergs/s, which im 



cluster halos. Also shown is the isothermal profile sometimes 
used to describe the cluster. In the following we adopt value 
CAe = 5, which is consistent with both observational con- 



plies that temperature values used by Pierpaoli et al. (2001 ) 
are on average 5% higher than those used by Finoguenov 



et al. (2001). Second, 13 clusters are contained in both sam- 



straints (Allen, Schmidt 



Fabian 2001 ) and theor etical 
predictions (Navarro et al. 1997, Bullock et al. 2001). We 
will ignore the small expected variation of ca^ with cluster 
mass. Based on figure hi we estimate the error induced by 
this choice to be of order of 10% if M500 is used to deter- 
mine M2oon„ • This error rises to 20-30% if M2500 is used 
instead to determine M2oon„ . 



3 RESULTS 

The observed M — T relation for T > 3kcV clusters relevant 



for abundance studies as given by Finoguenov et al. (2001) 
is 



M: 



2oon„ 



= 10^ 



'/i-'(l ±O.2)M0 ( - 

V/ 



T 



(5) 



pies. Linear regression with zero offset gives again the same 
5% deviation. The two methods give consistent results and 
we apply this small correction in the analysis below. An- 
other small effect which goes in the same direction is to cor- 
rec t back the redshift depe ndence of T which was applied 
by [Finoguenov et al. (200ll ) in the M - T fitting. Together 
these corrections make an effective (3 — 1.9 for the fiducial 
model. 

We emphasize that T is used here just as a label and 
does not necessarily need to be the correctly calibrated tem- 
perature of the cluster, as long as the mass assigned to this 
T is correct and there is a monotonic relation between the 
two. For this method to be unbiased we must only address 
the bias in the mass determination. However, the masses 
are derived from the assumption of hydrostatic equlibrium 
combined with the measurements of the temperature and 
gas profile. For a polytropic profile T oc p^as^ ^^id beta- 



where /3 = 1.75±0.25 and we used the relation between M500 
and M2oonm based on the adopted cluster profile (figure |]|). 
The value of /3 quoted above is for the fiducial model with 
Sim ~ 0.35, but in fact it varies little with ilm, increasing 
only by 5% for Qrn ~ 1- This is because while Ac in equation 
^ increases with f2m the mass M2oonm extrapolated from 
M500 decreases, compensating each other for a profile that 
is close to isothermal. 

While the scaling in equation ^ is in agreement with 
the theoretical predictions the normalization of the M — T 
relation is not. If we used Ma^ as done in previous work we 



model for X-ray intensity profile Jx oc 
one finds 

..3 

M(r) oc r/3x7 



k rc ) 



l/2-3/3x 



(6) 



Here Vc is the core radius of the beta-model profile. We see 
that the derived mass is linear in /3x, 7 and T (note that 
the /3x here should not be confused with (3 as a parametriza- 
tion of M — T relation), each of which has some uncer- 
tainty associated with it. Uncertainty in mass determina- 
tion from surface brightness profile fitting is expected to be 
small. Somewhat more controversial are T profiles derived 
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from ASCA, which have bee n interpreted by s ome groups 
(e.g. Finoguenov et al. 2001, Markevitch 1998) to show a 



decline in T at larger radii, in agreement wi th theoretica l 
predictions,while others do not find this (e.g. White 200C ). 
Such d eclin es are not seen in general from Chandra (Alien 
et al. 2 )01), although the radii probed are smaller and this 
is not unexpected. Bepp oSAX does see a decline in most 



of the clusters observed (De Grandi & Molendi 2001) and 



for the few clusters tha t overlap with the ASCA sample in 
Finoguenov et al. (2001 ) the drop in T out to rsoo is in a rea- 
sonable agreement between the two. It is also in agreement 
with XMM-Newton for the clusters that have been observed 
by bot h satellites, althoug h the XMM-Newton sample is still 
small ( Arnaud et al. 2001 ) . The measure d average polytropic 
index 7 1.2 ( |Finogu cnov et al. 2001) is close to theoret- 
ical predictions (Komatsu & Seljak 2001) and we verified 



that t his also provides a reasona ble fit to the BeppoSAX 
data (De Crandi & Molendi 2001) in the outer parts of the 



cluster where most of the mass is (in the inner parts the 
polytropic assumption fails, since it cannot predict a sharp 
flattening of the temperature profile at a radius larger than 
the core radius, as seen with BeppoSAX data). 

To derive the constraints on the power spectrum am- 



pUtude I we compute tne cumulative 1 function at an ettec- 



tive mdaii redshifl of ihe sample s: 0.05. We compare it 



to match the best fitted model of Pierpaoli et al. (2001) 
around T = 6.5keV, which is the pivot point where most 
of the sensitivity is. The results are not very sensitive to 
this choice since the shapes of cumulative mass function are 



very similar. Our simplified approach, which avoids a proper 



likelihood evaiution, is reasonable as long as the estimated 
errors in the two methods are comparable. Large errors can 
introduce bias since for a steeply declining T function low 
T clusters can scatter into high T and produce a cumula- 
tive T f unction w ith a higher amplitude in this exponential 
regime ( Pen 1995 ) . The scatter in Al — T relation is compa- 
ra ble or perhaps somewh at larger than the estimated errors 
in Pierpaoli et al. (2001). In the latter case this would lead 



to a slight overestimate of erg. Here we will assume the errors 
are comparable and we do not correct for this effect, which 
would further reduce erg. 

Figure ^ shows the comparison between the best fitted 
theoretical prediction for f2m — 0.35, Qa = 0.65, F = 0.2 



model and the data sample of Pierpaoli et al. (2001). Here 



r is the shape parameter of the power spectrum. The best 
fitted value for this model is as = 0.7. A more general fit 
around this model is erg = 0.7(n,„/0.35)~'''"' (F/0.2)''°*, 
which is supposed to give a few percent accuracy over 
0.2 < < I and 0.1 < F < 0.3. Note that the dependence 
on Q,m is less steep than the usual Q^'^, a consequence of 
the different mass temperature relation and different mass 
function and mass definition. We estimate the error by vary- 
ing /3 by one standard variation in both directions. This gives 
erg = 0.7±0.06 for the fiducial model. We emphasize that this 
is only the statistical error. The systematical error, which is 
much harder to estimate, is discussed below. 



4 DISCUSSION AND CONCLUSIONS 

The value for erg ~ 0.7 ± 0.06 for = 0.35 is about 25% 
lower than the most recent value derived from the local clus- 



10-i 



.,10-' 



10-' 




10-= 



Figure 2. Cumu 



kT(keV) 



Pier- 



I ative temperature function constructed by 
paoli et al. (2001) compared to the best fitted model for a flat 
cosmology with Qrn = 0.35. 



ter sample (Pierpaoli et al. 2001). While it is very clo se to the 



value obta ined from high redshift cluster catalog ( Borgani 



et al. 2001 ), it is important to note that they use f3 = 1.2 in 
their analysis and so if the actual value is /3 ~ 2 as suggested 
by the local M500 — T observations then their value for erg 
would be even lower. For the same value of /3 the discrep- 
ancy in the normalization between the local and the high 
redshift sample is at a more than 3-er level and difficult to 
explain as a statistical fiuctuation. This points to perhaps 
additional systematic effects beyond the M — T normaliza- 
tion discussed here that may be affecting either the local 
or the high redshift sample. We also note that our lower 
value for erg agrees well with other recent attempts to cir- 
cumvent the use of simulat ions by adopting the measured 
masses from weak lensing ( Viana, Nichol, fc Liddle 2001 ) 
and hydrostatic equilibrium ( iReiprich fc Boehringer 200 1[ ). 



The latter analysis in particular is very close in spirit to our 
approach. 

Even though the method used here does not rely on sim- 
ulations and so in principle should be more reliable, it still 
contains several possible systematic uncertainties. We men- 
tioned some before: the mass determination from spatially 
resolved temperature maps has perhaps a 20% uncertainty. 
It is difficult, although perhaps not impossible, to invoke 
systematic errors in the mass determination to bring the 
observed M — T relation in agreement with the adiabatic 
simulations results. We know however that additional pro- 
cesses must be included in simulations. Cooling, which is 
relatively well understood, although numerically challeng- 
ing, makes a significant difference in the normalization of 
M — T relation and its inclusion brings the predicted re - 
lation closer to the observations (Muanwong et al. 2001). 
Underestimation of errors in M — T can also bias the re- 
sults and if these are larger than assumed here then erg will 
be reduced further, although probably by not more than a 
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few percent unless the scatter in M — T relation is seriously 
underestimated. 

On the theoretical side the extrapolation from M500 to 
has roughly a 10% uncertainty in m ass. The use of the 
mass function by 3hcth fc Tormen (199i^ ) with Ma^ rather 
than Jenkins et al. (2001) mass function with M^oosi^ would 
decrease the value of ag by 7%. Rather than to construct 
mass function based on a spherical overdensity at a virial 
radius, which is not directly observable, it is more useful to 
focus on mass functions using masses defined within some 
smaller radius, such as rsoo, where mass can be directly mea- 
su red. A step in this d irection has been undertaken recently 
by Evrard et al. (2001), who evaluate mass function at A'hoo, 
which is only 40% above M500 for ca^ = 5. Unfortunately 
they do not provide a universal mass function valid for all 
models, but a fitting formula for the two cosmological mod- 
els they simulate. We find that for fiat Qrn = 0.3 mode l their 



mass fu nction a grees with ours if the mass definition of Jenk 
ins et a . (2001) is used together with ca„ = 5. Similarly, M. 



White (private communication) has evaluated mass func- 
tions using M500 for several models and again we find the 
agreement with our mass function is very good. This can be 
further improved in the future by defining the mass as close 
as possible to the observational definition (e.g. in projection 
for the case of lensing and X-rays) . The agreement between 
different groups gives us confidence th at the mass function 
has converged to the form proposed by Jenkins et al. (2001) 
and this should not be a major source of uncertainty in the 
future. On the other hand, use of different mass functions 
can explain up to 10% differences in derived as in the past 
work. 

Our results are actually good news for CDM models, 
which are facing significant problems on small scales (see 



Alam et al. 2001 for an overview). The models predict too 
steep density profiles of halos and too many subhalos within 
halos, both of which are related to predicting too much small 
scale power. If the lOMpc normalization can be significantly 
reduced this not only reduces the amplitude of the power 
spectrum, but also allows for a stronger tilt, which can fur- 
ther reduce the power on scales below IMpc where most of 
the problems are. Such tilted AC DM model was recently in- 
vestigated by Alam et al. (2001). It was shown there that 
the halo structure problems with CDM are significantly alle- 
viated and the tilted low as model gives halo concentrations 
in line with observations. Subhalo abundances are reduced 
as well, but not to the point where there would be too few 
to mach the observations once photoionization suppression 
is included ( [Bullock 20m\) . 

Perhaps the most important message of the present pa- 
per is that the systematic errors associated with the cluster 
abundance normalizati on are sti ll significantly larger than 
often assumed (see also Voit 2000 ) and so one should be cau- 
tious when ruling out cosmological models based on this con- 
straint. Our results suggest that even models with as ~ 0.6 
3 are not excluded at this point. Although at 
present the cluster normalization may be less reliable than 
previously thought, the prospects for the future are more 
promising. Instead of concentrating on calibration of M — T 
relation from the simulations, which may be uncertain for 
some time in the future, one should focus on obtaining better 
data to empirically calibrate the relation through the mass 
determinations from lensing, velocity dispersions and appli- 



cations of hydrostatic equlibrium. The latter should be espe- 
cially promising as well calibrated spatially resolved temper- 
ature profiles become available with Chandra, BeppoSAX 
and XMM-Newton. These mass measurements should be 
tied to the temperature determinations of clusters used in 
the construction of cluster sample and compared to mass 
functions constructed using masses which are observation- 
ally defined. As larger samples become available one will 
be able to construct mass selected samples directly with- 
out the need of using temperatures at all. Ultimately this 
approach should provide a reliable constraint on the cosmo- 
logical models at lOMpc scale. 
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